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Abstract: Hydrological models are often linked with other models in cognate sciences to understand the 
interactions among climate, earth, water, ecosystem, and human society. This paper presents the 
development and implementation of a decision support system (DSS) that links the outputs of hydrological 
models with real-time decision making on social-economic assessments and land use management. 
Discharge and glacier geometry changes were simulated with hydrological model, water availability in semi- 
arid environments. Irrigation and ecological water were simulated by a new commercial software MIKE 
HYDRO. Groundwater was simulated by MODFLOW. All the outputs of theses hydrological models were 
taken as inputs into the DSS in three types of links: regression equations, stationary data inputs, or dynamic 
data inputs as the models running parallel in the simulation periods. The DSS integrates the hydrological 
data, geographic data, social and economic statistical data, and establishes the relationships with equations, 
conditional statements and fuzzy logics. The programming is realized in C++. The DSS has four remarkable 
features: (1) editable land use maps to assist decision-making; (2) conjunctive use of surface and 
groundwater resources; (3) interactions among water, earth, ecosystem, and humans; and (4) links with 
hydrological models. The overall goal of the DSS is to combine the outputs of scientific models, knowledge 
of experts, and perspectives of stakeholders, into a computer-based system, which allows sustainability 
impact assessment within regional planning; and to understand ecosystem services and integrate them into 
land and water management. 


Keywords: decision support system; hydrological modeling; ecosystem services; land management; socio-economic 
indicator; Tarim River Basin 


1 Introduction 


The debate over the effectiveness of integrated water resources management (IWRM) in practice 
has lasted for years (Jeffrey and Gearey, 2006; Biswas, 2008; Quevauviller, 2010; Giordano and 
Shah, 2014). As the complexity and scope of IWRM increases, the difficulties of hydrological 
modeling are shifting from the model itself to the links with other cognate sciences and to 
understand the interactions among water, earth, ecosystem and humans. This paper presents the 
development and implementation of a decision support system (DSS) that links the outputs of 
hydrological models with real-time decision making on social-economic assessments and land use 
management. The DSS assists the decision-making in a qualitative manner based on the outputs of 
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hydrological models and knowledge of experts in cross-disciplinary fields. 

There are currently many DSS models with economic assessments (Wenkel et al., 2013; West 
and Turner, 2014; Yue et al., 2014; Pedro-Monzonis et al., 2016), focusing on single topic (e.g., 
climate change, land use, agriculture, water account), and many socio-economic models (Flérke et 
al., 2013; Arnell and Lloyd-Hughes, 2014; Herrero et al., 2014; Visconti et al., 2015) without 
simulation of hydrological processes. Few model is capable of integrating meteorological, 
geographical, ecological, social and economic factors, which is based on the simulation of 
hydrological models. Little research has been done to show the interactions of so many cross- 
disciplinary researches on TWRM to understand ecosystem services (ESS) and integrate them into 
land and water management. Due to model complexity issues, many DSS models have to consider 
less issues within a specific research framework (McCown, 2002; Basso and Ritchie, 2015). 
However, the model complexity of DSS is worth the effort (Chenoweth et al., 2004; Power et al., 
2015) to increase model accuracy on the comprehensive analysis in IWRM. 

Since the last decade the Chinese government has been promoting the development of western 
China. The demographic development and socio-economic change has led to a rapid change of land 
use systems in the Tarim River Basin, where it is the habitat for more than 80% of Populus 
euphratica trees in China, and has substantially affected the quantity and quality of arable soil, 
surface water, and groundwater (Zhao et al., 2009; Zhang et al., 2014; Rumbaur et al., 2015). These 
changes in soil and water have large impacts on the crop production and natural vegetation (Chen 
et al., 2015; Zhang et al., 2015), such as the P. euphratica trees. 

Since 2011, the German Ministry of Science and Education Bundesministeriums fuer Bildung 
und Forschung (BMBF) established the Sino-German SuMaRiO (www.sumario.de) project for 
sustainable management of river oases along the Tarim River. A cross-disciplinary consortium of 
11 German and 9 Chinese universities and research institutes joint together for the research on 
SuMaRiO and the DSS. Project SuMaRiO focus on realizable management strategies, considering 
social, economic and ecological criteria. This will have positive effects for nearly 10x10° 
inhabitants of different ethnic groups. The DSS is the main outcome of SuMaRiO. The overall goal 
of the DSS is to integrate all crucial research results of SuMaRiO, including stakeholder 
perspectives, into a model based decision support system, which allows a sustainability impact 
assessment within regional planning on land and water management in the Tarim River Basin to 
understand ESS and interactions with water, earth and humans. The DSS is an indicator based tool 
that enables stakeholders and decision-makers to evaluate the consequences of their intended 
actions, which helps to implement sustainable land and water management measures in the 
upcoming development plans (Disse, 2016). 


2 Materials and methods 


2.1 Study area 


This research is conducted in the Tarim River Basin. With the mainstream of 1321 km and located 
at the northern edge of the Taklimakan Desert, the Tarim River is the longest inland river in China. 
The Tarim Basin is known worldwide for its natural resources, extreme arid climate and vulnerable 
ecosystem. In this region, agricultural water consumption and allocation management are crucial 
to address the conflicts among irrigation water users from upstream to downstream (Yu et al., 2015), 
which cause severe water scarcity problems on ecosystem (Liu and Chen, 2006) and sustainable 
land management (Feng et al., 2001; Liu and Meng, 2011). 

This research focused on the mainstream of Tarim River, starting from Aral till Taitema Lake. 
With a mean annual precipitation of around 40 mm and more than 3000 mm potential 
evapotranspiration per year, the region is extremely arid with vulnerable ecosystem. Water in the 
mainstream of Tarim River is mainly provided by high mountain glacier melting, which feed the 
Tarim River through its tributaries. Water resources and ecosystem stability are the most obvious 
and sensitive issues in this hyper-arid region. 

The oases along the Tarim River is covered by riparian forest, grassland, reed, shrubs, and 
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farmlands. The main crop of the irrigated fields is cotton, which is also a major income source for 
the local people. The Tarim River Basin has a global reputation of producing high quality cotton, 
which provides more than half of the cotton production in China. 


2.2 Research fields 


The interdisciplinary research challenges are clustered in five interrelated work blocks, WB 1-WB 
5 (Fig. 1), with WB 1 on organization with stakeholders and management issues, WB 2 on regional 
climate change and discharge of the Tarim River tributaries, WB 3 on sustainable water and land 
use management, WB 4 on ecosystem services (ESS) and ecosystem functions (ESF), and WB 5 
on multi-level socio-economic assessment of ecosystem services and implementation tools. The 
DSS will identify realizable management strategies, considering social, economic and ecological 
criteria. 


Global assessment of land use dynamics on greenhouse 
gas emissions and ESS 


WB 5: Multi-level socio-economic 
assessment of ecosystem 
services and implementation tools 


WP 2.1: Monitoring and modeling SuM; aRiO 5.1: Multi-level economic system assessment 


of cryosphere B, 
WP 2.2: > a : 5.2: Transdisciplinary assessment of ESS for 
2.2: Regional climate scenarios and urban areas regarding dust and heat stres: 
ium-term forecast of precipitation 
WP 2.3: Climate change impact 
on water discharge 


5.3: Actor-based decision support 
for land and water management 


WB 1: Organization 

WP 1.1: Project coordination and 
equipment management 

WP 1.2: Scenario management 


WP 1.3: Stakeholder dialogue and 
coordination of knowledge transfer 


WP 1.4: GIS and DATA management 


WP 4.1: ESS and ecosystem functions of riparian 


MEST Lilly ener eres pei ecosystems and aquatic biodiversity 
P 3.2: Hydrology, salinity and biomass WP 4.2: ESS and ccosystem functions of 
production on the local scale (10 km?) non-irrigated land use systems 
P 3.3: Upscaling to the regional scale (200 km?) Wp 4.3: ESS and ecosystem functions of urban 
P balance and peri-urban oasis vegetation 
00 km 


Fig. 1 Work blocks with interdisciplinary research fields (Disse, 2016). WB, work block; WP, work package; 
ESS, ecosystem services. 


Research content of the DSS covers the following aspects in the Tarim River Basin: water 
resources, climate, biodiversity, demographics, energy consumption and production, poverty and 
health, economic development, land management, and scenario management. A number of these 
research fields involve more than one WB, and thus need expert knowledge from interdisciplinary 
cooperation. 


2.3 Indicators of DSS 


Under scenario assumptions, possible actions and their impacts are estimated in a quantitative way 
with the help of sustainable indicators (Table 1). Climate indicators, socio-economic indicators and 
management indicators are input indicators, and ESS indicators are output indicators. Input 
indicators give users opportunities to change the scenario on the baseline or perspective in the 
planning years, and output indicators demonstrate the simulation results caused by the management 
alternatives. 

User can define the inputs of DSS according to the reference values, possible ranges and certain 
rules. Simulation years can be chosen from 2012 to 2050. The climate projections were based on 
three emission scenarios, nine global climate models, and additional two regional climate models. 
Four types of climate scenarios, CCLM RCP 2.6, CCLM RCP 4.5, CCLM RCP 8.5, and REMO 
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A1B, are included in the DSS. For socio-economic indicators and management indicators, users 
can change the reference values on the management alternatives in the planning years. The results 


of their actions will be demonstrated by the output indicators. 


Table 1 


Climate indicators 


Socio-economic indicators 


Input and output indicators of the decision support system (DSS) 


Management indicators 


Temperature rising (°C) 


Precipitation increase (%) 


Cotton production costs (CNY/hm?) 


Average costs for the production of 


fruits (CNY/hm?) 


Drip irrigation share in total agricultural 
irrigation area (%) 


Subsidy for drip irrigation (CNY/hm?) 


Input Average costs for the production of 6 $ TETE 
mda other crops (CNY/hm?) Household and industry water allocation (%) 
Cotton price (CNY/t) Flooding of natural vegetation (%) 
Average price of fruits (CNY/t) 
Average price of other crops (CNY/t) 
Ecosystem ESS ESS indicators 
Provisioning services Cotton production (x106 t) 
Provisioning services Fruit production (x10° t) 
Agriculture 
Provisioning services Production of other crops (x10° t) 
Provisioning services Farmer's income (x10° CNY) 
Provisioning services Biomass production (x106 t) 
Regulating services Drifting dust control by riparian forest (kg) 
Output Regulating services Sand mobilization control by riparian forest (x10° t) 
ee Riparian Forest 
indicators 


Grassland 


Regulating services 
Regulating services 
Supporting services 
Provisioning services 
Provisioning services 
Regulating services 


Regulating services 


Wind control (Attenuation in % at 2 m height) 
Carbon storage (x106 t) 

Mean species (x10° t) 

Apocynum production (x10° t) 

Reed production (t) 

Drifting dust control by grassland (kg) 


Sand mobilization control by grassland (x10° t) 


Note: Apocynum is used as a medicinal plant and a fiber crop in the arid region of Central Asia. ESS, ecosystem services. 


2.4 Logics of DSS 


Equations, conditional statements, and fuzzy logic consist the three logical relationships, which 
connect the parameters in the DSS. 


2.4.1 Equations 


Equations are formed by expert knowledge, theorem, empirical equations and literature findings to 
build the relationships among the parameters in the DSS. In sum, there were 115 equations 
established in the DSS. Two examples are given in the following text. 

Due to the absence of industry in the region, water consumption in the catchment is comprised 
of domestic, flooding of natural vegetation, irrigation water use and water losses (Eq. 1). Inflow of 
the first sub-catchment is given by a model of water availability in semi-arid environments (WASA) 
depending on the climate scenarios. Because no bifurcation exists in the mainstream of Tarim River, 
inflows into the sub-catchments are calculated in a simple way each month. 


Q(n +1) = O(n) — O(n) x P(domestic) — O(n) x P(flood) x k — IW x 0.964 , (1) 


where Q(n+1) is the inflow into the next sub-catchment (m°), O(n) is the inflow of the current sub- 
catchment (m°), P(domestic) is the interpolated value of domestic water use (%), P(flood) is the 
interpolated value of flooding of natural vegetation (%), k is the flood distribution over the month 
(%), and IW is the irrigation water demand of current sub-catchment (m3), which multiplied with a 
loss factor of 0.964. 

Cotton production is calculated cell-by-cell, then aggregated in sum for the DSS output (Eq. 2). 
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ws 1-61-22) x Ymxm/100, (2) 
ET 


c 


where Y. is the cotton production (t) in each cell, ET, is derived from Penman-Monteith method 
and influenced by crop factors, ET, is calculated by MIKE HYDRO, which identified the actual 
rate of crop evapotranspiration under the effects of soil water stress, Ym is the maximum harvest 
yield (kg), m is the cotton production distribution over months (%), and Ky is the yield response 
factor representing the effect of reduction in evapotranspiration on lost yield. Ky values were 
obtained from FAO Irrigation and Drainage Paper No. 33 (Doorenbos and Kassam, 1979). 


2.4.2 Conditional statements 


The knowledge of experts are also used to formulate the 'if-then' rules to solve modeling problems. 
Such conditional construct perform different computations whether a condition is evaluated to be 
true or false. Conditional statement is a convenient method if the rules are clear and homogeneous. 
Three examples are given as follows. 

For a grassland cell, if groundwater level < —5 m in 9 months each year and continues for 7 years, 
then land use type of the cell will become unused land in the 8" year. 

For a high density riparian forest cell, if groundwater level < —10 m in 9 months each year and 
continues for 7 years, then the grid cell will become riparian forest with low density in the 8" year. 

For a low density riparian forest cell, if groundwater level < —10 m in 9 months each year and 
continues for 7 years, then the grid cell will become unused land in the 8" year. 


2.4.3 Fuzzy logic 


In the DSS, some relationships cannot be precisely expressed as equations or conditional statements 
and a useful method under such circumstance is called fuzzy logic. Fuzzy logic is a logical approach 
to deal with uncertainty management (Zadeh, 1983). Elements have degrees of membership in the 
fuzzy set. Membership functions represent the degrees of truth in vaguely defined sets. Unlike 
possibilities or conditions, the truth values of variables may be any real number between 0 and 1. 
Examples of tree height, crown area and drifting dust control by riparian forest are demonstrated 
in Tables 2-4. 


2.5 Models linked with DSS 
Discharge and glacier geometry changes were simulated with hydrological model WASA (Giintner 


Table 2 Fuzzy logic of tree height 


Membership function for tree height (m) 


Status Left bound Middle Right bound 
Low 3 = 5 
Medium 5 7 12 
High 12 - 20 
Fuzzy rules for tree height 
Rules Groundwater level (m) Flooding (m) Tree height (m) 
If High And High Then High 
If Low And Low Then Low 
If Medium And Medium Then Medium 
If High And Medium Then High 
If Medium And High Then Medium 
If Low And High Then Medium 
If High And Low Then Medium 
Note: "—" means not necessary. 


Table 3 Fuzzy logic of crown area 
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Membership function for crown area (m°) 


Status Left bound Middle Right bound 
Low 2 — 6 
Medium 5 7 12 
High 10 = 20 
Fuzzy rules for crown area (m°) 
Rules Groundwater level (m) Flooding (m°) 
If High And High Then High 
If Low And Low Then Low 
If Medium And Medium Then Medium 
If High And Medium Then High 
If Medium And High Then Medium 
If Low And High Then Medium 
If High And Low Then Medium 


Table 4 Fuzzy logic of drifting dust control by high density and low density riparian forest 


Membership function for drifting dust control by riparian forest (kg/hm”) 


Status Left bound Middle Right bound 

Low 0 = 3.3x10° 
Medium 3.3x10° 4.2x10° 5.8x10° 

High 6.7x10°° - 8.3x10° 


Fuzzy rules for drifting dust control by riparian forest high density 


Rules Gronn aa Tree ere € Drifting dust control by riparian forest 


high density (kg/hm*) 
If High And High Then High 
If Low And Low Then Medium 
If Medium And Medium Then High 
If Low And Medium Then High 
If Medium And Low Then High 


Fuzzy rules for drifting dust control by riparian forest low density 


Drifting dust control by riparian forest low 


Rules Crown area (m7) Tree height (m) eid ohn) 
If High And High Then Low 
If Low And Low Then Low 
If Medium And Medium Then Low 
If High And Medium Then Low 
If Medium And High Then Low 


and Bronstert, 2004). The calibration and evaluation of the model were conducted in the headwater 
catchments of Tarim River. WASA provides daily discharge input into the DSS. Based on different 
temperature rise and precipitation change projections, four climate scenarios, A1B, PCP2.6, 
RCP4.5, and RCP8.5, could be selected by users (Duethmann et al., 2016). 

For the surface flow, hydrological model MIKE HYDRO (Yu et al., 2015) provides inputs for 
the DSS on a catchment level. MIKE HYDRO is a comprehensive deterministic and physically 
based modeling tool for the simulation of water flow, water supply/demand, soil moisture and crop 
growing. It has an integrated modular structure with basic computational modules for hydrology 
and hydrodynamics. The DSS and MIKE HYDRO have the same delineation of sub-catchments 
(SC): SC1 from Aral to Xingiman, SC2 from Xinqiman to Yingbazar, SC3 from Yingbazar to Qiala, 
and SC4 from Qiala to Taitema Lake. The irrigation water demand, ecological water for the natural 
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vegetation, seepage losses, fruit production, and crop yields were calculated and summarized as 
inputs into the DSS. 

For the groundwater simulation, a MODFLOW model was computed with internal links to the 
DSS on a 500 mx500 m cell level (Fig. 2). MODFLOW is running parallel with the DSS in the 
simulation. The daily water head of each cell is simulated in MODFLOW, then updates the cells in 
the DSS. Groundwater recharge is considered via four different sources: river leakage, irrigation 
water seepage, infiltration during flood season, and ecological water percolation in the lower 
reaches. 


N \ J 


Input 
discharge 


=> 


scl 
r4 


/ Surface flow 


Groundwater recharge 
boundary conditions 


rs \ Groundwater 
x Linked on cell level |\OQDFLO 


500 mx500 m per cell 


500 mx500 m per cell 


Fig. 2 Hydrological models linked with decision support system (DSS). Model of water availability in semi-arid 
environments (WASA) provides discharge. MIKE HYDRO provides water consumption and water balance. 
MODFLOW runs parallel with the DSS on groundwater simulations. SC, sub-catchment; SC1, sub-catchment from 
Aral to Xingiman; SC2, sub-catchment from Xingqiman to Yingbazar; SC3, sub-catchment from Yingbazar to Qiala; 
SC4, sub-catchment from Qiala to Taitema Lake. 


It is a challenge to find links between the DSS and other models. The key issue is that the outputs 
of different models are not uniform as for the input into the DSS. The parameter sets and their 
relationships vary largely from model to model. The outputs of different models can either serve as 
stationary data inputs, regression equations, or dynamic data inputs into the DSS as the models 
running parallel in the simulation periods. 

As MIKE HYDRO linked with the DSS on the sub-catchment level, data and equations are 
bridging the two models. MIKE HYDRO provides simulation results of ecological water for the 
natural vegetation in the upper and middle reaches in the DSS. Crop productions are linked by yield 
equations calculated from evapotranspiration. Seepage losses were calculated in MIKE HYDRO 
and aggregated in monthly values into the DSS. Irrigation water consumptions were distributed in 
MIKE HYDRO and integrated on sub-catchment level in the DSS. 

Particularly, MODFLOW has internal links with the DSS. As the same delineation on the cell 
level, both models are running parallel in the simulation periods. At the beginning of each year, 
water heads from MODFLOW results are updating the groundwater levels in the DSS. 
Groundwater plays a crucial role in the ecosystem, and has large influence on a number of output 
indicators in the DSS. 


2.6 Graphical user interface of DSS 


A user-friendly graphical user interface (GUI) was developed to assist the decision-makers and 
common users with Chinese and English versions available at the moment. Labels and instructions 
allow users interact with the DSS more conveniently. The input indicators have default values, 
acquired from database, literature and expert knowledge. Calculation years start from 2012 to 2050. 
The land use map is editable as well in the baseline year, to assist decision-making on land use 
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changes. The GUI is designed to be user-friendly so that more stakeholders and scientists could get 
involved in using and improving the model. Implementation of the DSS bridges the gap between 
research and practice through interviews, workshops, and feedbacks. 


3 Results and discussion 


3.1 Land use changes 


Users can freely define the input land use types by clicking the cells in the GUI. In this trial, land 
use types are remained as they were in the default map which acquired from MODIS data. The 
baseline year is 2012, and all the values of input indicators are kept business as usual. the planning 
years are 2030 and 2050, which indicate land use changes in the near and far future (Fig. 3). 


er i eS ee ee (ere ee ae 


a a 


Legend E Riparian forest (high density) [iy eee 
E hy - 
CO Unused land E Riparian forest (low density) a 
E Tarim River Urban vegetation ee a 
E Grassland E Other agriculture S, 
2050 C © Natural vegetation ™ Cotton ere 
| E Natural vegetation otto 3 


Fig. 3 Land use changes from 2012 to 2030 and then to 2050 


Grassland, natural vegetation and forest will suffer varying degrees of deterioration in the near 
and far future if the scenario is kept business as usual. In the upper, middle, and lower reaches, the 
decay of green areas is obvious. Similar results were shown by the latest findings (Liu et al., 2016; 
Ly et al., 2016). Because of the low precipitation and flow rate in the dry seasons, vegetation is 
dependent on summer flood and groundwater. Additionally, riparian forests would be gradually 
vanished in some regions by 2050, especially in the middle reaches. Since cotton and other 
agricultural fields are determined by anthropogenic activities, farmland areas remain unchanged 
throughout the simulation period. However, the consequences of farmland changes can be revealed 
by output indicators if a user changes the land use types in the first place. In the middle reaches, 
where riparian forest (mostly P. euphratica) is mainly growing, forest degradation will be severe 
by 2050 (Fig. 4). 

Simulation results indicate that both forest area and density have been reduced in various degrees. 
In several hotspots, riparian forest areas have substantially dwindled or even faded away. Theses 
hotspots are all located on the north of the river. Population is one reason to identify a hotspot, 
where deforestation could be a big threat to the living conditions of local residents. Another reason 
is the vanishing points of forest in several regions, which would be a destruction to ecological 
balance in arid land. Groundwater recession is the reason for forest decay in the model, as 
deforestation has been prohibited in the oases along the Tarim River. Continuous low groundwater 
level for more than 7 years is considered fatal to the forest. The correlations between groundwater 
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and P. euphratica were also discussed by Thomas (Thomas et al., 2017), who made further 
discussions on the threats for these precios trees. 


aak ANTS E: pals 
Hotspots of forest degradation ' 
Fig. 4 Forest degradation in the middle reaches of Xinqiman to Qiala sub-catchments from 2012 to 2050 


3.2 Downstream outflow 


The natural terminal of Tarim River is Taitema Lake. However, due to long-term water interception 
downstream, the outflow of Tarim River is studied at Daxihaizi Reservoir. The reservoir is located 
358 km upstream of the river till Taitema Lake. The outflow of Tarim River in future years is 
dependent on the upstream discharge and water consumption along the oases, which were 
calculated from the other hydrological models. Simulation results indicate an increasing trend of 
outflow in 2020, and decreasing trends in 2030 and 2040 (Fig. 5). 


1500 
—-~— 202 —< 2020 —@®— 2030 —e— 2040 


1000 


500 


Outflow («10° m3/month) 


-500 
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 


Fig.5 Downstream outflow in the Daxihaizi Reservoir 


Peak volumes of outflow appeared in July, with 7.3x108, 9.5x108, 6.9x108, and 5.9x108 m? in 
2012, 2020, 2030, and 2040, respectvely. From Figure 5, it was noticeable that summer flood last 
two months in July and August, and drop dramatically in September. These results agreed with the 
outputs of hydrological model (Yu et al., 2015). Dry season occupied most of the year. In June, 
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water deficits occurred in all the simulation years. Particularly, near 4.0x108 m? of water shortage 
was predicted in June 2040. A noticeable water deficit also showed up in September 2040, which 
is the first time that the outflow water is on shortage in September. Upstream reservoirs need to 
prepare enough storage for the water crisis in future. 


3.3 Socio-economic outputs 


Cotton production, fruit production, biomass production, mean species of plants, farmer's income, 
wind control, drifting dust control by riparian forest, drifting dust control by grassland, and sand 
mobilization control by grassland are shown as examples for the socio-economic assessments (Fig. 
6). These indicators were calculated separately in the four sub-catchments of Aral to Xinqiman, 
Xingiman to Yingbazar, Yingbazar to Qiala, and Qiala to Taitema Lake. Simulation period is from 
2012 to 2050, and all the inputs were kept as default values. 
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Cotton production has an increasing trend in all the sub-catchments. The largest cotton 
production region is within SC1, but the fastest production growth comes from SC4. This result 
complies with the study of Xu et al. (2008), who warned that the agricultural development in the 
lower reaches may cause water scarcity problems for the ecosystem. Fruit production also has a 
rising trend, and develops rapidly in the lower reaches. Biomass production experiences an 
increasing period till 2018, then drops in 2019 and remains steady afterwards. The reason for this 
sudden change may be because the rules in DSS are not complex and smooth enough, and 
quantitative accumulations lead to a qualitative transformation in 2019. 

The number of mean species of plants has a steady growth, similar to farmer's income. Because 
of the development of economic crops, farmer's income in the lower reaches has the largest growth 
among all the sub-catchments. Due to lacking of clear policy and data in future years, change of 
population and migration are not considered in current version of the DSS. Wind control and 
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drifting dust control by riparian forest remain stable, even the forest area in shrinking. This is 
caused by the contribution of growing tree height and crown area. In comparison, drifting dust and 
sand mobilization control by grassland have declined in varying degrees until 2050. 

One scenario was tested as an example for land use management and assessment (Fig. 7). In the 
baseline year of 2012, 4x100 cells (1x10* hm?) of cotton fields are changed into grassland in the 
first sub-catchment to investigate the consequences of this action in future years until 2030. 

The simulation results indicate the reduction of cotton production and farmer's income, and 
increase of sand mobilization control by grassland after land use change. Cotton production 
provides the main income for local farmers. From a local decision-maker's point of view, if the 
losses of farmer's income are too large to bear, then this land use change scenario may not be 
accepted. However, if sand mobilization control is considered more important than the economic 
losses, the decision-maker may adopt this scenario, and mitigate the farmer's losses in other ways. 
Additionally, the increased portion of sand control is shrinking in the future years. This 
phenomenon indicates either groundwater level is too low in the fields, or groundwater salinity is 
too high for grass to grow. 
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Fig.7 Editable land use map in the DSS. In SC1, 4x100 cells of cotton fields are changed into grassland in 2012. 
Cotton production, farmer's income and sand mobilization control by grassland are investigated as output indicators 
of this action. 


3.4 Uncertainties and further developments 


In the DSS, the delineation of regions is relatively homogeneous in multiple criteria, including 
hydrological, meteorological, geographical, ecological, socio-economical, and political aspects. 
Therefore, the complexity of the management practices cannot be completely reproduced. 
Equations, conditional statements and fuzzy logic all have limitations in the delineation of 
parameter relationships, and thus cause uncertainties by temporal and special changes on specific 
matters. Socio-economic and management indicators (e.g., cotton price, farmland area and drip 
irrigation coverage) are often sensitive to the change of government policies and water authority 
regulations, which could cause large uncertainties in future planning years. 

Although the hydrological models were calibrated and validated to increase the possibility of 
more reliable results, the DSS has not been calibrated so far. However, a comprehensive model 
calibration requires the feedback from local stakeholders, more statistical data, and a better 
understanding of the DSS itself. Since the DSS involves knowledge on varying aspects, no 
individual developer is able to master all the criteria. Interdisciplinary cooperation is as important 
in the implementation and optimization phase as it was in the system developing phase. 

It is difficult to say which scenario has the most robustness, or one scenario is more preferable 
than the other. Based on a number of trials and comparisons with real data, most outputs of the DSS 
were reasonable and comply with the knowledge and perceptions of local stakeholders. However, 
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the sensitivity of soil salinization change are lower than expected. Because winter irrigation uses 
extra water to washes down the salt, salinization on the surface soil should be changed obviously 
in winter. The process of winter irrigation is not included in the simulation of DSS at the moment, 
and can be considered as a special stress to increase system reliability in future improvements. 


4 Conclusions 


A model based DSS was developed to assist the stakeholders with decision-making on sustainable 
land management and socio-economic assessments in the oases along the Tarim River. The 
development and implementation of DSS involves knowledge of experts in interdisciplinary 
research aspects, as well as the experiences and feedbacks from a large number of local 
stakeholders. The DSS has notable features in a number of research aspects. 

(1) Editable land use maps in the GUI were developed to assist stakeholders of decision-making 
on land use management. The clicking and drawing on the land use map provide stakeholders a 
catchy approach for examining land management alternatives. If the scenario is kept business as 
usual in the simulation, then natural vegetation, riparian forest and grassland would suffer varying 
degrees of deterioration in the near and far future. In the upper, middle and lower reaches, the decay 
of nature vegetation are obvious. 

(2) The conjunctive use of surface and subsurface water resources provide more accuracy on the 
system. The DSS integrates the hydrological data, geographic data, social and economic statistical 
data, and establishes the relationships with equations, conditional statements and fuzzy logics. 
Under scenario assumptions, possible actions and their impacts are estimated in a semi-quantitative 
way with the help of climate indicators, socio-economic indicators, management indicators, and 
ESS indicators. Socio-economic outputs illustrate more crop productions and farmer's income, but 
weaker sand mobilization and drifting dust control in future. 

(3) The integration of expert knowledge on interdisciplinary studies gave specific insights on the 
interactions among water, earth, ecosystem and human activities. The system does not only 
consider direct impacts, but also side effects among different matters. For instance, increased 
irrigation water may raise crop production, but it can also cause less water for winter flooding, 
therefore aggravate salinization on the field, and lead to less yield in the end. 

(4) Links with hydrological models help reduce system complexity and increase model accuracy 
of the DSS. MODFLOW is running parallel in the simulation periods with the DSS. WASA 
provides discharge inputs, and MIKE HYDRO simulates irrigation and ecological water 
consumptions. All the outputs of theses hydrological models were taken as inputs into the DSS. 
After the calculations of water balance in the DSS, simulation results indicate an increasing trend 
of outflow in 2020, and a decreasing trend in 2030 and 2040. 

Though many expert knowledge and stakeholder feedbacks were considered in the developing 
phase, the DSS should not be recognized as an end product. Lots of estimations and uncertainties 
still remain in the system, which requires further research and development. The GUI was designed 
to be user-friendly so that more common stakeholders can use and improve the system, and make 
the steps in bridging the gap between research and IWRM practice. 
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